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Abstract. We have studied numerically the random interchange model and related loop 
models on the three-dimensional cubic lattice. We have determined the transition time 
for the occurrence of long loops. The joint distribution of the lengths of long loops is 
Poisson-Dirichlet with parameter I or 


1. Introduction 

The random interchange model is a stochastic process where transpositions are selected 
at random. The product of these transpositions gives a random permutation and the main 
question deals with its cycle structure. We consider variants where possible transpositions 
are restricted to nearest-neighbours of a regular cubic lattice. We also consider related loop 
models where “crosses” are replaced by “double bars”. We provide evidence that a phase 
transition takes place where macroscopic loops occur. We give good estimates of the values 
of the parameters at the transition point, which should help to better comprehend the model. 
Finally, we compute moments of the lengths of the loops; it turns out that they are identical 
to those of the Poisson-Dirichlet distribution. 

The random interchange model was invented by Harris El; Toth used it as a represen¬ 
tation of the spin ^ quantum Heisenberg ferromagnet m- Angel [3] and Hammond m 
obtained results for the model on trees. Schramm considered the variant on the complete 
graph and he proved that the joint distribution of the large cycle lengths is Poisson-Dirichlet 
of parameter 1, following a conjecture of Aldous [15]; see also [4] for some simplifications and 
extensions. Alon and Kozma have obtained remarkable identities that give the probability 
of cyclic permutations in terms of eigenvalues of the graph Laplacian, for arbitrary graphs 
[2]. These have allowed Berestycki and Kozma to obtain further results on the complete 
graph [5|. 

A similar model with “double bars” instead of “crosses” was introduced by Aizenman 
and Nachtergaele in order to describe the spin | quantum Heisenberg model and the spin 
1 model with biquadratic interactions [T]. It should be noticed that the representation of 
quantum systems involves the extra factor with ^ = 2,3... In the survey [8|, the 

authors conjectured that the joint distribution of the lengths of long loops, in dimensions 
three and higher, is given by the Poisson-Dirichlet distribution of parameter 0 . The two 
representations of Toth and Aizenman-Nachtergaele were recently combined so as to describe 
quantum models that interpolate between the two Heisenberg models, such as the spin ^ 
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quantum XY model and further spin 1 models with SU(2)-invariant interactions [18]. For 
these models, the joint distribution of long loops should be Poisson-Dirichlet with parameter 
^/2. The consequences of this structure have yet to be worked out. One such consequence 
is to identify the nature of symmetry breaking in the spin 1 model m- 

The Poisson-Dirichlet distribution is conjectured to be a common feature of “loop soups” 
in dimensions three and higher. This was confirmed numerically in lattice permutations |9] 
and in 0(N) loop models [13]. This was also confirmed, with a mathematically rigorous 
proof, in an annealed model of spatial permutations [6]. A numerical study of the spin 1 
model 120 ] also provides indirect evidence, as explained in m- However, this conjecture is 
far from being accepted nowadays. It thus seems necessary to verify it also in those loop 
models that are related to quantum spin systems. 

We introduce the random loop models in Section |2.1[ The conjectures about the joint 
distribution of the lengths of long loops are explained in Section [2^ Numerical evidence is 
presented in Section [^ for the occurrence of a phase transition; we also discuss a comparison 
with bond percolation and quantum Heisenberg models. We investigate the presence of the 
Poisson-Dirichlet distribution in Section |4j 

2. Random loop models 

2.1. Definitions. Let A = {1,..., C and let Sa denote the set of edges (nearest- 
neighbours) in A. Let ^ > 0 and u G [0,1]. To each edge of Ea is associated an independent 
Poisson point process on the interval [0,/3] with two kinds of events: 

• crosses occur with intensity u; 

• double bars occur with intensity 1 — u. 




Figure 1. Graphs and realisations of Poisson point processes, and their 
loops. In both cases, there are exactly two loops. 

See Fig.[^for an illustration. Given a realisation cj, we define loops by moving in the vertical 
direction, and by jumping to the neighbour whenever a cross or a double bar is encountered. 
If it is a cross, one continues in the same vertical direction; if it is a double bar, one continues 
in the opposite vertical direction. We assume periodic boundary conditions in the vertical 
direction. All trajectories must close, resulting in loops. We use the following notation for 
the relevant random variables. 

• I/ 2 (ci;),... denote the vertical lengths of all the loops in decreasing order 
(repeated with multiplicities); notice that 0 < Lj{uj) < ;d|A| for all j. 

• ^i(co’), • • • denote the “shadow lengths” of the loops in decreasing order (re¬ 

peated with multiplicities); that is, G {1, 2,..., |A|} is the number of sites at time 
0 in the jth loop. 
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Notice that for all realisations cj, we have Lj{uj) = /3\A\ and = 1^1- The 

following are random partitions of [0,1]: 

( Lijuj) L2H \ W \ ,,,, 

V m ’ ^iAi V |A| ’ |A| ^ ^ 

When (3 is small, crosses and double bars are scarce and loops are small. But a phase 
transition occurs as P grows and some loops become macroscopic. A few random partitions 
are displayed in Fig. they were measured in a cube of volume 160^, for a value of {3 that 
is above the critical parameter. Occurrence of macroscopic loops is manifest. 


H - 1 - 1 I I iii^ 

0.00 0.25 0.50 0.75 1.00 



0.00 0.25 0.50 0.75 1.00 



0.00 0.25 0.50 0.75 1.00 


Figure 2. Samples of random partitions observed in a cube of size L = 160, 
for R = 1 and /3 = 1. Notice that m{P) appears to be constant with 
approximate value 0.9. 


2.2. Phase transition and universal behaviour. The patterns of the random partitions 
suggest a very interesting strong law of large numbers. It is worth describing in details, 
since it is expected to occur in all models of “loop soups”. 

Conjecture 1. There exists m{P) sueh that, as L ^ oo, we have for almost all realisations 


K 


lim lim 


K —>-CXD L —>-CXD _ ^ I I 


hm hm \ = 1 — mi^). 

k^ooL^oo ^ A 
j:£j<k ' ' 


( 2 . 2 ) 


The function m{/3) represents the mass of points in long loops. It is equal to 0 when (3 is 
small, and it becomes positive when (3 crosses the transition point. Eq. (2.2) says that loops 


are either microscopic (the length is of order 1) or macroscopic (the length is of order |A|); 
the number of sites that belong to loops of intermediate, mesoscopic length, has vanishing 
density. 

Conjecture 1 is also expected to hold with Lj/p instead of with the same m{P). 
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The next conjecture is about the joint distribution of the lengths of the macroscopic 
loops, which should be Poisson-Dirichlet (PD). Let us recall the closely related Griffiths- 
Engen-McCloskey (GEM) distribution and its “stick breaking” construction. Let Xi, X 2 ,... 
be i.i.d. Beta(l,'d) random variables; their probability density function is 'd(l — for 

0 < s < 1. The following is a random sequence with GEM('d) distribution: 

(Xi, (1 - Xi)X2, (1 - Xi)(l - X2)X3, ...). (2.3) 

It is not hard to verify that the sum of all these numbers is 1 with probability 1. Rearranging 
the numbers in decreasing order, we get a random partition with PD[op]('d) distribution. 
Multiplying each element by m, the distribution is PD[o,rn] ('^)- We can formulate the con¬ 
jecture about the joint distribution of macroscopic loops. 


Conjecture 2. For any k, the joint distribution 0 /^ 1 ,... converges to the joint distri¬ 
bution of the first k elements of a random partition with PZ)[o, 7 n(/ 3 )] W distribution^ where 




ifu = 0 or 1; 
ifu G (0,1). 


Here, m{fi) is the same quantity as in Gonjecturej^ In random loop models with weights 
^#ioops^ the conjecture holds with'd = if r = 0,1, and d = 0/2 \i u ^ (O 5 !)• 

The heuristics for this conjecture goes back to Aldous’ ideas for the complete graph, which 
Schramm eventually managed to turn into a proof m- Its relevance for models with spatial 
structure was suggested in [HI [9] . In summary, the idea is to consider the stochastic process 
restricted on random partitions. Adding crosses or double bars result in splits or merges of 
elements of the partition. The rate at which two long loops merge may seem at first sight 
to depend on the exact geometry of the loops, which is very intricate. But averages take 
place when the loops are macroscopic, which is the case in dimension 3 and higher. The 
effective stochastic process is then a standard split-merge process, whose invariant measure 
is Poisson-Dirichlet. 

The addition of a transition (cross or double bar) between different loops always results in 
a merge. When u = 1, or when r = 0 on a bipartite graph, the addition of a transition within 
the same loop always results in a split. It follows that the Poisson-Dirichlet distribution has 
parameter'd = 1. When 0 < r < 1, the addition of a transition within the same loop may 
split it, or it may rewire it. The splits occur therefore at half the rate of the merges, and 
the Poisson-Dirichlet distribution has parameter ij = ^. See [HI El [El H] for more detailed 
explanations. 


3. Phase transition and critical parameters 


3.1. Fraction of sites in long loops. We seek a convenient expression for the mass of sites 
in macroscopic loops The expressions of Gonjecture turn out to be inconvenient 

and we use Gonjectn re instead. It allows to relate m{/3) with the moments of the lengths 
of the loops, see Eq. ( |3.4[ ) below. We now give a derivation of this result. 

We start with 


E 




loop 

j>l x,yEA 


1 


loop , 


x,y£A 


(3.1) 
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If we accept Conjecture the probability that two sites that are far apart, belong to 
the same loop, is given by the probability that they both belong to long loops — this is equal 
to — times the probability that two random numbers in [0,1] belong to the same 

element in the partition. The latter probability can be calculated using the GEM random 
sequence (2.3). The probability that both random numbers belong to the jth element is 

/ dSi / dS2 IEgEM('i?) (lsi,S 2 eith element) = 1 EgeM('i?)(^^) 

Jo Jo 

= IEBeta(l,i?) ((1 - X)‘^y IEBeta(l,i?) {X‘^) . 

Elementary computations give EBeta(i,i?) ((1 “ and EBeta(i, 7 ?) = (, 9 +i)%+ 2 ) • 

The probability that two distant sites belong to the same loop is therefore approximately 
equal to 

2 _ m{Pf 


(3.2) 


P((a:,0) O (2/,0)) =TO(^)2y](-^^)^ - 


i>i 


(i9 + l)(t? + 2) id + 1 


(3.3) 


Using (13, we obtain an expression that is convenient for numerical calculations, namely 


j>l 


(3.4) 


Our numerical results are displayed in Eig. As expected, m{P) is zero for small [3 and 
positive for P large enough. m{P) is also continuous and increasing. Notice also that, when 
R = 1, it converges to 1 as /3 ^ oo; indeed, all sites belong to long loops. It converges to a 
value smaller than 1 when r = 0 or | because a density of small loops remains present in 
the system. 


3.2. Value of the critical parameter /3c( r). The previous results confirm the existence 
of a phase where m(/3) is positive. We define the critical parameter by 

/3c(r) = inf{/3 : m(/3) > 0}. (3.5) 

It is instructive to estimate its value and to investigate its dependence on the parameter u. 
Our numerical results are depicted in Eig. We find that I3c{u) is a convex function of u. 
Its minimal value occurs when u is close to 0.5. The derivative of Pc{u) diverges at r = 0 
and u = 1. In retrospect, this is perhaps not surprising: The model with weight has 

more symmetry at r = 0 and r = 1, SU(2) rather than U(l), so that a minor change in the 
value of u has major consequences. This also explains why Pc{\) < /3c(0) and /3c(|) < /dc(l): 
Less symmetry means reduced fluctuations so that spontaneous magnetisation occurs more 
easily in systems with U(l) symmetry rather than SU(2). Somehow, this explanation should 
retain validity when the weight is not present. 


3.3. Comparison with bond percolation, Ising, and quantum Heisenberg models. 

There are natural comparisons between the critical parameters of the model of random loops, 
of the Ising model (a.k.a. the q = 2 random cluster model), and the quantum ferromagnetic 
and antiferromagnetic Heisenberg models. 

A. Bond percolation. Given a realisation uj of crosses and double bars on x [0, there 
corresponds a percolation configuration r] = r]{uj) where rjxy = 1 if at least one transition 
occurs on xy x [0, /3], and rjxy = 0 otherwise. The percolation parameter is p = ^{rjxy = 1) = 
1 — e~^ . Let C{uj) denote the set of percolation clusters of One can check that each 
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Figure 3. Numerical values for (y^)^) as function of /3. The size 

of the cube is L = 80 in the top row and L = 160 in the bottom row. The 
parameter u takes values 0 (left), 0.5 (centre), and 1 (right). The points 
represent the values of many realisations; the curve gives the 

average. 


loop 7 G C{(jj) is contained in some cluster C G C{uj) in the sense that 7 C C x [0,/d]. Then 
|C(u;)| < \jC{uj)\ and the presence of long loops is possible only when percolation occurs. 

The critical parameter for the cubic lattice is Pc = 0.2488, which gives = — log(l — 
Pc) = 0.286. We have indeed found that ^c{u) > /d^®^ for all u G [0,1], and the inequality is 
strict. 


B. Random cluster and Ising model. The random cluster model is similar to bond 
percolation, but with configurations 77 receiving the extra weight The random 

cluster model is closely related to the q state Potts model; their transition parameters 
satisfy = 1 — g = 2 is equivalent to the Ising model with 

2 ^ising _ ^Potts(g)^ Section 3.8] for an excellent introduction to this topic. The Curie 

temperature for the three-dimensional Ising model is Tc = 4.51, which allows to deduce that 




,r.c.(g=2) 




Potts (g=2) 


= 0.443. 


There is no natural comparison between (31'^' of the random cluster m odel , and f3c{u) 


(2) 

of the random loop model. But it can be compared with /d^®^ and Pc 7 see (3.7). 


C. Quantum Heisenberg models. Let 5'^, 5'^, 5'^ denote the usual spin operators in 
that satisfy [5'^, 5'^] = id's and further cyclic relations. Consider the quantum Hamiltonian 
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u 


Figure 4. Value of the critical parameter Pc as function of u. 


that acts on the Hilbert space (8)a:eAC^: 

Ha = -2 Y. {SISI + {2U-1)SISYSX). (3.6) 

{x,y}cA 
\\x-y\\ = l 

The case u = 1 corresponds to the spin | Heisenberg ferromagnet; the case u = ^ is the 
quantum XY model; the case r = 0 is unitarily equivalent to the Heisenberg antiferromagnet, 
provided A is bipartite. 

The partition function and the quantum correlations can be expressed using random loop 
models where realisations uj receive the extra weight [T6l[Tl|T8]. In particular, P^\u) 

is equal to the inverse Curie temperature of the quantum model. Numerical studies have 
found that = 0.59 [IT] and = 0.53 [14]. 

The extra weight encourages the system to have more, smaller loops, so we can expect 
Pc {u) > Pc{u) for all u. Besides, we observe that Pc TO) < Pc Tl), as in the absence of 
the weight. 

Regarding the quantum XY model, Stefan Wessel has just performed numerical calcu¬ 
lations and he obtained the value Tc = 1.008 ±0.001 for the Hamiltonian with interaction 
~^{x,y}(^x^y ^x^y) EB* This implies that Pc^\^) ~ 0.496. We conjecture that P^\u) 

has a shape similar to that of Pc{u). 

Let us summarise the discussion above with the following inequalities; for all u G [0,1], 


'Pc{u) G [0.313,0.361]' 

/^r.c. (9=2) 


/3P®'' = 0.286 < 


= 0.443 


< e [0.496,0.59]. 


(3.7) 
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4. Joint distribution of the lengths of long loops 


4.1. Calculation of the moments of Poisson-Dirichlet. We check the presence of the 
Poisson-Dirichlet distribution by looking at its moments, following m- Let Til > ''' > Uk 
be integers. The calculation of the moments can be achieved by starting from another 
representation of the Poisson-Dirichlet distribution due to Kingman [12]. Let Zi,..., Zjv 
be i.i.d. random variables with Gamma(^) distribution (that is, their probability density 


/r(^) for 0 < s < oo). Let S = Zi • 


+ Ztv- Consider the sequence 

(4.1) 


f Zi Zn \ 

VS'’*"’ S J 

and reorder it in decreasing order, so it forms a random partition of [0,1]. As N ^ oo, this 
partition turns out to converge to PD[op]('d). The following two observations are keys to 
our calculations: 

• S' is a Gamma('d) random variable; 

• S is independent of (^,..., ^). 

For given integers ni,..., >0, using the independence of S from the partition, we have 

N\ rZk\'^^ 


E 


'PD 


[ 0 , 1 ] 


(i?) 


( s 


yni 


distinct 


■ A J {N -k)l 


lim 


Euff 


= lim 


A^! 


(f) ) 


(4.2) 


N^oo{N-k)\ E(S'"i+-+"fe) 

J. A^! r(i?)E(Zp ...zp) 
N—^oo (AT — k^\ “1“ fii “1“ ■ ■ * “1“ 

We also used E(S'“) = r(?? + a)/r(??). Since the ZjS are independent, 


E(zp...zp)=n 

Recall that r{'d/N) ~ N/'d as N 


Ti'd/N + m) 


E 


PD[o,il(-?) 


T{^/N) ■ 

00 , so that (jv-fc)r(,?/jv)fc We obtain 

i9'=r(i9)r(ni)...r(nfe) 




11 vdfc>i 
distinct 


r('d + ni H-h n/c) 


(4.3) 


(4.4) 


This important formula appears in [13]. Its derivation there is different; it involves another 
loop soup model, assumes the presence of Poisson-Dirichlet, and uses a “supersymmetry” 
method. 

Eq. (4.4) holds for Poisson-Dirichlet on the interval [0,1]. The formula for the interval 
[0,m] is identical, except for the additional factor . Combining Conjectureand 

Eq. ( |4.4[ ), we get the following exact formula. 

Conjecture 3. The moments of the lengths of the loops are given by 


E 




jiv,ifc>i 

distinct 


#r(^)r(ni)...r(nfe) 

r('d + Ri + • • • + R/c) 


where 'd = 1 for u = 0 or u = 1, and 'd = ^ for 0 < u < 1. 

With minor modifications, this formula also applies to a wide range of “loop soup” models. 
And indeed, it was derived in [13] in the context of 0(N) loop models. 
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4.2. Numerical results. We now calculate numerically some moments of the joint distri¬ 
bution of the loops and compare the results with Conjecture First, moments involving a 
single loop. Let be the value of m(/3) in Conjecture l^when ni G N and rij = 0 for 

j >2. Explicitly, we have 




r(^ + m) 

'dr{'d)r{ni) 



(4.5) 


Our numerical results deal with /3 = 1, u E {0, |,1}, and ni G {2,3,4,5} and they are 
listed in Table They show that mnj^{(3) is quite constant in rii, apart from numerical 
fluctuations and finite-size corrections. This confirms Conjecture]^ Notice that it confirms 
in particular that i) is either 1 or ^, depending on the value of u. 



0 

u 

1 

2 

1 

2 

0.8925 

0.9585 

0.9310 

3 

0.8968 

0.9587 

0.9276 

ni ^ 

0.8815 

0.9595 

0.9217 

5 

0.8930 

0.9528 

0.9356 


Table 1. Numerical values of for j3 = 1, and L = 160. 


Next, let be as mmiP), but with n 2 7 ^ 0. The expression is 


mn,,n2{P) 


r(i9 -I- ni -1- ng) 
t?2r(i9)r(ni)r(n2) 




l/(ni+n2) 


(4.6) 


The numerical values are listed in Tables [2||^ for r = 0, 1. Notice that we avoid the values 

= 1 because of undesirable effects due to small loops. 



2 

n2 

3 4 

5 

2 

0.9031 

0.8924 0.8915 

0.9001 

3 


0.9056 0.9001 

0.8985 

ni ^ 


0.8987 

0.8947 

5 



0.8946 


Table 2. Numerical values of mn^^ri 2 {P) for r = 0, /d = 1, and L = 160. 
Here, mL{P) = 0.872. 


We observe that mni,n 2 (/^) does not depend much of ni, n 2 , as expected from Conjecture 
This also guarantees that the value of'd has been conjectured correctly. Variations in the 
values of 771 ^ 1,712 (/d) can be dismissed as random fluctuations and finite-size effects. 

To summarise, we have studied numerically the joint distribution of the lengths of macro¬ 
scopic loops in a family of loop models in three dimensions. These loop models are motivated 
by their close relations to certain quantum spin systems and include in particular the ran¬ 
dom interchange model. We have observed the presence of the Poisson-Dirichlet distribution 
with parameters'd = 1 and ^ as conjectured in mm- 
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2 

n 2 

3 4 

5 

2 

0.9576 

0.9579 0.9562 

0.9593 

3 


0.9671 0.9600 

0.9580 

ni ^ 


0.9539 

0.9479 

5 



0.9523 


Table 3. Numerical values of = 0-^? /^ = 1? L = 160. 

Here, mL{P) = 0.949. 



2 

n 2 

3 4 

5 

2 

0.9286 

0.9295 0.9318 

0.9315 

3 


0.9350 0.9256 

0.9341 

ni ^ 


0.9349 

0.9282 

5 



0.9210 


Table 4. Numerical values of mni,n 2 (/^) for r = 1, /3 = 1, and L = 160. 
Here, mL{P) = 0.925. 


Acknowledgments: We are grateful to Stefan Wessel for sending us the result of his 

numerical calculation of the critical temperature of the quantum XY model. 
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